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^3 . ABSTRACT 



We present an improvement of global Metropolis updating steps, the instanton hits, used in a 
hybrid Monte Carlo simulation of the two-flavor Schwinger model with staggered fermions. These 
hits are designed to change the topological sector of the gauge field. In order to match these hits 
to an unquenched simulation with pseudofermions, the approximate zero mode structure of the 
J> ' lattice Dirac operator has to be considered explicitly. 

o 

oo . 

O ■ Introduction 

, The Schwinger model^ (SM) shows a lot of features of great interest, such as topological sectors of 

"j^ I the gauge field and corresponding zero modes of the Dirac operator according to an index theorem.^ 
Closely related to that is the anomalous breaking of the chiral C/(l)-symmetry.^ With these properties 
Ph! the SM may be used as a toy model for QCD, or with a coupled scalar field, for the electroweak 
(— I I theory."^ 

For a study of the mentioned non-perturbative features in those theories lattice calculations with 
- Monte Carlo (MC) methods are needed among others. However, while topological field configurations 
r> ' are studied because of their effect on the fermionic sector, the reverse relation, i.e. the influence of 
\ dynamical fermions on the dynamics of topological quantities, is usually not treated within direct 
lattice calculations. One reason is the large computational effort for systems with dynamical fermions 
close enough to the continuum limit. 

In Ref. 5 we studied the topological properties of the lattice SM with staggered fermions using 
the hybrid Monte Carlo (HMC) algorithm with pseudofermions.^ As staggered fermions describe two 
flavors of Dirac fermions in the continuum limit, we compared the results with continuum calculations 
of the two-flavor SM on the two-dimensional torus 72^'^ We found good qualitative reproduction even 
of the role of the topological sectors and zero mode contributions. However, we took a rather small 
lattice of size 6 x 16. One reason was that the HMC algorithm, which is often used for calculations 
with dynamical fermions, does not tunnel the potential barriers between different topological sectors 
for large /?. The additional instanton hits,^ designed to solve this problem in pure QED2, do not 
work on larger lattices. 

We present here a matching of this instanton hits to dynamical fermions treated with the pseudo- 
fermion method. In particular, whenever the topology of the gauge field is changed, the zeromode 
structure of the Dirac operator changes simultaneously. This could be taken into account with local- 
ized instanton hits, which are optimized in the pseudofermion background. The improved instanton 
hits overcome the bad large volume behavior. 
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Before we come to this improvement, we give a brief description of the topological properties of 
the lattice SM and its implications for a MC simulation with pseudofermions. Some useful results of 
Ref. 7 on the zero mode structure in the continuum SM are given in the appendix. 

1. Topological Properties of the Lattice Schwinger Model 

For the lattice formulation of the SM with staggered fermions we start with the action (x = 



S[U,4>A]=Sg + Sf 



xy •ry , 



x,y 



^'^P^l,x U[x,n]^x+e^„y f^[jy^^]<^a;'-e^,y 



(1.1) 

(1.2) 



Px is the Wilson plaquette, /? is 1/e^ in lattice units. For the fermion field we took periodic boundary 
conditions.^ Since the SM is superrenormalizable, the continuum limit a ^ is achieved for large 
P DC 1/a^. In Ref. 5 we found that for /? = 10 the continuum results are reproduced with rather small 
lattice effects. The transition to the strong coupling behavior is found at /? ~ 1. 

The lattice definition of the topological charge Qtop can be derived from the lattice field strength 
Fx byio 

1 



Fx e [- 



■vr,7r 



Q 



top 



(1.3) 



This charge has integer values and coincides with the Chern index k in the continuum limit. Corre- 
sponding to the index theorem Eq. ( A. 3 ) in the continuum, one finds 2|(5fop| approximate zero modes 
(AZMs) for large enough (3. Exact zero modes arise only in gauge field configurations of zero measure 
in configuration space. Therefore we may write down the naive formula for fermionic four-point 
functions 



< 



>-- 
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V[U] e 



-Sg[U] 



det(M[?7]) 



(1.4) 



In the continuum version of such expectation values on a finite volume the zero modes must be 
treated explicitly, leading to zero mode contributions from the topological sectors k = ±1,® see also 
Ref. 12. It was one of the aims of Ref. 5 to identify such contributions in lattice simulations. 

There are two difficulties in this task, which shouldn't be mixed up. Firstly we expect low values 
of the fermion determinant det(M[f7]) in the topological sectors Qtop = k, due to AZMs for large 
p. Therefore the weights qk of those sectors are small. On the other hand the AZMs should lead to 
large values of certain combinations of the fermion propagator leading to contributions of the sectors 
Qtop = il) which do not vanish in the continuum limit. In order to gain sufficient statistics in theses 
sectors, it is convenient to use a modified action, in which the sectors Qtop pushed. The 
weights Qk can then be determined from the modified weights q'j^,^ see also Ref. 13. 

The second difficulty arises from the potential barriers between the different topological sectors. 
The HMC algorithm did not manage to tunnel these barriers for large values of p. In Ref. 5 we used 
the so-called instanton hits,^ working well in pure QED2. They are given as Metropolis updating 
steps with a certain global proposal. However, with dynamical fermions (i.e. with pseudofermions) 
the acceptance rate for these proposals breaks down rapidly for lattices larger than the 6 x 16-lattice 
taken there. 



^ It turns out that antiperiodic boundary conditions for the fermion fields can be absorbed in a shift of the constant 
part of the gauge field. Thus the choice of fermionic boundary conditions doesn't matter in fluctuating gauge fields. 
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In order to improve these instanton hits one has to consider the imphcations of the zero mode 
structure in an algorithm with pseudofermions. In the pseudofermion method the weight factor for 
the gauge field configurations is substituted by 

det(M[;7]) = J P[$^«>] e-^f[^'*l, Sp[U,^] = <^>\M[U]M^ [U])-'^^ . (1.5) 

Here the pseudofermion field $ is a complex bosonic auxiliary field. For staggered fermions it can be 
taken only on the even sites of the lattice, thus avoiding the additional spectrum doubling due to the 
product M[?7]M''^[C/].^^ So we must only consider the even AZMs of this product, i.e. we have IQtopI 
even AZMs x'iU] for large /?. The HMC algorithm now generates configurations [U,^] determined 
by the action S = Sg + Sp in Eqs. (pTiP, ( |1.5| ). As this action contains M[U] to a negative power, the 
AZMs lead to a strong suppression of configurations with <I> not orthogonal on [U] 

\<^,X'[U]>\^0 . (1.6) 

2. Instanton Hits with Pseudofermions 

The instanton hits proposed in Ref. 9 are metropolis update steps of the gauge field with a 
special global proposal, designed to change the topological charge Qtop by ±1. Detailed balance of 
such Metropolis steps can be most easily attained by a symmetric proposal probability Pp 

Pp{U, $ ^ U', = Pp{U', ^' U, ^) (2.1) 

and an acceptance probability for these proposals in the Metropolis decision 

Pa{U, $ ^ U', $') = rmn{l, e"^^) , AS = Sg[U'] + Sp[U', - Sg[U'] - Sp[U, $] . (2.2) 

The proposal of [U' , <!>'] is implemented by an interpolating gauge field AC/ with Qtop = 1 

U^p] = AC/[,,^] • C/[,,^] , K = '^.- (2.3) 

The proposal probability Pp is symmetric, if AU does not depend on U, and if the inverse change 
{AU)~^ is proposed with equal probability. The computational cost of these instanton hits is es- 
sentially given by the determination of S[U',^'] for the Metropolis decision Eq. ( |2.2| ), including a 
Conjugate Gradient (CG) solution of the linear equation {M[U]M^ [U]) X = ^. The success of these 
hits is limited by the mean acceptance probability Pa- 

In pure QED2 the interpolating gauge fields are determined by an optimization of P4. Since we 
want to choose AU independent of U, we assume the most likely starting configuration U, i.e. the 
minimum of the old sector Qtop = n. The aim is to propose U' such that the action is minimal in 
the sector Qtop = nit 1. Thus AU interpolates between two neighboring minima, it is an (euclidean) 
instanton. The minima are given by a constant field strength 

= 27rQtop/{LiL2) , so AF, = ±27r/(LiL2) . (2.4) 

With dynamical fermions the situation changes completely. Since the pseudofermions are held 
fix, the gauge of the interpolating field AU is no longer arbitrary. The choice of a minimal gauge 
for AU is not unique, because the necessary kink in a topologically non-trivial configuration can 
be situated everywhere on the lattice, see e.g. Ref. 15. As well, Sp[U,^] depends strongly on the 
constant part of U, the toron field tp, see Eq. ( |A.2| ). In Ref. 5 it was therefore necessary to combine 
the instanton hits with simultaneous toron shifts 

^x,/.] ^ e^^*" > At^ = R27T/L^, (2.5) 
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Table 2.1: Mean acceptance rate Pa and mean reduced acceptance rate Pr for Qtop = — >■ ±1 
proposals on a 12 x 12-lattice, /? = 4. 





old (global) instanton hits 


local instanton hits (r = 6) with x-opt. 


Pa: 


0.0005(1) 


0.0043(3) 


pR-- 


0.0262(8) 


0.1040(14) 



where i? is a random number between —1/2 and 1/2. However, also with this slight modification, 
the acceptance rate Pa breaks down rapidly for larger volumes. 

For large enough /? and thus sharp AZMs, the main part in AS" comes from the AZMs in the 
new gauge field configuration XnewW]- They are no longer orthogonal on the pseudofermion field 

leading to a high Sp[U'] and thus to low acceptance rates Pa in Eq. (^I^). In order to control 
this effect we measured the acceptance rate Pa, and the reduced acceptance rate Pr for proposals 
Qtop = ^ Qtop = if- In the definition of the reduced acceptance rate, the effect of the lowest 
mode is separated 

Pr ^ mzn(l, e-^^«[^'*l) , Sn[U, f ] = S[U, $] - | < <1>, xo > l' Aq"^ , (2.6) 

where xo is the lowest normalized even eigenmode of M[C/]Af^[C/], with eigenvalue Aq- It appears 
that the strongest suppression of the acceptance rate is in fact due to the zero mode part of ASp, 



see Table 2.1. In this sense we want to achieve 

<'^,X^oid>-0 ^ <<i>,xiew>~0 (2.7) 



for as many as possible new AZMs xiiew^ see Eq. (|2.6|). Let us therefore consider the zero mode 
structure in more detail. 

In the continuum SM there is a multiplicative relation between the even zero modes in the gauge 
field backgrounds >l^(x), AA^{x), and A'^{x) = A^{x) + AA^(x), with Chern indices n>0, 1 and 
n+1, respectively, see Eq. ( |A.8| ). If the zero mode x in A.A^{x) is essentially constant on most of the 
torus 

X^(x)~ca;^forx0^, (2.8) 
it follows from this multiplicative structure, that n of the n + 1 zero modes Xnew in are approx- 



imately given by the old zero modes x''oid everywhere, but in the region A. With Eq. ( [A.8D we can 
estimate 

<^,X'new>< C<^,X^M> + ^ < |$ | , | (x - C w) | > , (2.9) 

where m is an upper bound for \Xold{^)\- On the lattice, the first part of the rhs of Eq. ( |2.9D vanishes 
approximately in the simulation, see Eq. ( |2.6| ). The second part is small if the region of strong 
deviations A in Eq. (^) is small. 

Thus, on the lattice we look for interpolating fields with an AZM showing the behavior of Eq. (|2.8| ). 
The region of strong deviations A should be as small as possible. We take local instanton configura- 
tions with link variables Af/j^ = exp(iAA|'^ ^j) 



for < x^ < r 

(2.10) 

elsewhere 
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1-A{x) 
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• □ 




- □ 



1 . □ 




Figure 2.1: Absolute value A{x) (plotted as 1 — A{x)) and phase ip{x) of the AZM x^{x)/Q^{x) on 
a 32 X 32-lattice, Q^{x) = 0.0463(— 0.0424i) for xi,X2 = even (odd), x is only defined on the even 
sites, for simplicity the values on odd sites are given by the neighboring even sites. 



The field strength then is F^; = 27r/r^ for < < r , = elsewhere. The (even) AZM x^{x) 
in this configuration can be seen numerically to be 



(2.11) 



r{x) = fixwix) , 

where /(x) ~ 1, if x is far enough from the region of non- vanishing gauge field ^[x,;^], see Figure p.l| . 
Q,^{x) = a^{a^^) for xi,X2 = even(odd) corresponds to the differential form coj in the continuum for 
a^^ = —ia^, see Eq. (A. 7). This relation gets approximated on the lattice for increasing r. Other 
choices of local configurations (C^[x,/i] = 1 for ^fj. ^ [0)''')) l^ad to zero modes with the same principal 
behavior, yet with bigger deviations. 
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Figure 2.2: The mean acceptance rates PA{r) (full) and Pj^{r) (dashed). 



Plotting the mean acceptance rate PA{f) of such local instanton hits shows the effect of the 



localization, see Figure 2.2. The change in the pseudofermion action Sp favors small r, as indicated by 
the mean fermionic acceptance rate P^{r), P^ = min(l, exp(5p — Sp) ). This effect gets compensated 
by Sg for r < ropt , ropt depends on /?. 

Up to now we have optimized the interpolating gauge field AC/ independent of the actual configu- 
rations in the MC simulation. However, it is possible to choose /SXJ dependent on the pseudofermion 
$, which is held fix during the update by Af/. For fixed r it remains the freedom to shift Af/^ in 
Eq. ( p. 10 ) round the torus by x 



By this we minimize the lattice analogue of < [<I>|, [(x — ccj)! > in Eq. (|27 

X : ^ |'^(3;)| ■ \x^'{x — x) — n{x — x)\ = min 



(2.12) 



(2.13) 



n{x) is defined in Eq. (pTI] ) 
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Figure 2.3: Flux diagram for the improved instanton hits. 



In our MC algorithm we implemented the new instanton hits as shown in Figure |2.3| . The 
number Njh of instanton hits (combined with a preceding heatbath for the pseudofermion field) per 
Molecular Dynamics update is set equal to the number of Molecular Dynamics steps, in order to 
obtain approximately the same computational effort for both parts of the algorithm. Since the AZM 

must be calculated only once for a given lattice, the effort to determine x is significantly smaller 
than the CG solutions of the appearing linear equations. 



The use of matching x for given $ can be seen in Figure \2.4\ . It appears that there is no breakdown 
of Pa for large lattices. The acceptance rates are still smaller than 1%. However, the computational 
effort for a single instanton hit, essentially given by one CG solution, is low compared with the effort 
for the HMC updates. They need Nmd CG solutions, Nmd is the number of Molecular Dynamics 
steps {Nmd = 10 to 60 for lattice sizes 6 x 6 to 32 x 32). In Table 2A it can be seen, that there is 



also an improvement of the mean reduced acceptance rate Pr, in which the effect of the zero mode 



is separated, see Eq. (2.6). This means an improvement of the instanton hits also with respect to 
the higher modes. 

We measured the weights qn of the topological sectors Qtop = n^n = 1,2. For that we used the 
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Figure 2.4: The acceptance rate Pa in dependence on the lattice size L x L: For local hits with 
matched x (full line), with arbitrary x (dashed), for global hits (dotted). 
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Table 2.2: The weights qi and q2 vs. the lattice size L x L for /3 = 4. 



L 


6 


8 


12 


16 


24 


32 




0.0122(3) 


0.0295(6) 


0.073(3) 


0.128(7) 


0.171(12) 


0.239(10) 


Q2- 





0.00006(3) 


0.0014(4) 


0.0080(18) 


0.025(7) 


0.062(7) 



Table 2.3: The lowest mean eigenvalues Aj of {M[U]M'f[U]). 



12 X 12, 


/? = 4: 


Ai 


A2 


As 


A4 


As 


Ae 


Qtop — 







0.0198(4) 


0.0351(5) 


0.0627(4) 


0.0865(5) 


0.129(1) 


0.159(1) 


Qtop — 


±1 




0.0038(3) 


0.033(1) 


0.056(1) 


0.093(1) 


0.124(1) 


0.160(2) 


Qtop 


±2 




0.0018(3) 


0.0051(4) 


0.056(7) 


0.098(4) 


0.122(5) 


0.164(5) 


18 X 18, 


/? = 


9: 


Ai 


A2 


As 


A4 


As 


Ae 


Qtop — 







0.0115(1) 


0.0170(2) 


0.0321(2) 


0.0409(2) 


0.0663(3) 


0.0783(3) 


Qtop — 


±1 




0.00060(3) 


0.0190(5) 


0.0271(5) 


0.0488(6) 


0.0619(7) 


0.0824(8) 


Qtop — 


±2 




0.00055(5) 


0.0016(1) 


0.032(1) 


0.045(1) 


0.064(1) 


0.080(1) 


32 X 32, 


(3 = 


4: 


Ai 


A2 


As 


A4 


As 


Ae 


Qtop — 







0.0008(1) 


0.0022(2) 


0.0042(2) 


0.0064(2) 


0.0094(3) 


0.0122(3) 


Qtop — 


±1 




0.00052(6) 


0.0018(1) 


0.0036(2) 


0.0055(3) 


0.0084(4) 


0.0118(4) 


Qtop — 


±2 




0.0004(1) 


0.0013(2) 


0.0035(6) 


0.0061(7) 


0.009(1) 


0.011(1) 



equality g„ = due to a charge conjugation symmetry. The results for /? = 4 and different lattice 
sizes are given in Table 2^ . For large volumes the suppression of the non-trivial sectors gets weaker 
significantly. As reason we assume that on the lattice the suppression by the only approximate zero 
modes for a given value of /3 can be overcome by phase space factors for large volumes. The finite 
volume behavior, characterized in section 1, is then achieved again by increasing /? and fixing the 
volume in physical units. 

We want to conclude with the following remarks: Firstly, for the change from Qtop = n to 
Qtop = n + 1, there is one additional new zero mode, which is no subject of the matching of the 
interpolating gauge field AU and therefore in general is not (approximately) orthogonal on In 
particular, for Qtop[U] = 0, a zero mode Xoid determining the only zero mode Xnew in the background 
of AU ■ U exists only for special toron values = 0. In fact, for approximately vanishing toron 
values, the improvement of the instanton hits worked better than for the general case. On the other 
hand we found the acceptance rates for hits changing Qtop = ±1 ^ ±2 to be of the same order of 
magnitude as those for Qtop = — > ±1 on the larger lattices, where the higher sectors are not so 
strongly suppressed. 

Secondly, the arguments leading to the improved instanton hits are strictly based on the case of 
combining an instanton (and its zero mode) with another instanton, thereby increasing the overall 
topological charge. On small lattices we didn't find any evidence of approximate zero modes corre- 
sponding to a pair of local instanton and anti-instanton, which does not affect the topological charge. 
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The number of AZMs, determined by a significant gap in the eigenvalues A* of (M[?7]M^ [[/]), is given 
by Qtop, see Table 2.3. The situation for /3 = 4 on a 32 x 32-lattice is less clear. For us it is an 



open question, whether this is due to AZMs connected to instanton anti-instanton pairs or due to 
the ambiguity in the definition of the topological charge for the relatively low value /? = 4. Such 
an ambiguity becomes stronger on larger lattices. Note that a simple combination of instanton and 
anti-instanton solutions is not possible, due to the necessary spin flip in such combinations, see the 
appendix. Eigenvalues for other parameters and in the quenched case can be found in Ref. 16. 

Finally, a different picture would arise, if the zero modes of topologically non-trivial update 
configurations were localized, i.e. c = in Eq. Old, local zero modes would remain approxi- 

mately invariant under such updates. The update configurations should then be matched in order to 
minimize the scalar product of the pseudofermion field and the new zero mode. 
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Appendix 

The two-flavor SM in the continuum on the torus 72 was studied by S.I. Azakov, H. Joos. '''^ They 
partly use the geometric formulation with Dirac-Kahler (DK) fermions.^'' This means the fermions 
are described by inhomogenous differential forms (/)(x, H)dx^ , H=%, 1, 2, 12. In two dimensions these 
can be transformed into two flavors of Dirac fermions ipab{x)-, a = 1,2 is the spinor index, b = 1,2 is 
the flavor index. Since there is a systematic lattice restriction of DK fermions to staggered fermions, 
the geometric formulation is useful for a comparison of continuum and lattice model, in our case for 
the comparison of continuum zero modes with approximate zero modes on the lattice. In this brief 
description of the continuum zero mode structure we will use the Dirac language as far as possible. 

We consider the zero modes Xabix) of the Dirac operator in the gauge field A^{x) 



faa' 



[df, - ieA^ix)] Xa'bix) = . (A.l) 
The U{1) gauge field on T2 can be written 

2TTk 

Af,{x) = df,a{x) + e/a^6(x) + + —e^^x^, (A.2) 

d^a{x) is a pure gauge, €^dub{x) is a gauge field in Lorentz gauge, is the so-called toron field. The 
last term oc e^Xu is a representative of the Chern class with Chern index k ^ 7j. It is periodic up 
to (non-trivial) gauge transformations, and it has a constant field strength, i.e. it is a topologically 
non-trivial classical solution of the free field equations. For A; = 1 it is called an instanton. 
According to the Atiyah-Singer index theorem for two fiavors of fermions^^ 

ni — Ur = 2k , (A. 3) 

the Dirac operator in such a background field has ni = 2k left-handed zero modes for k > 0, rir = 2k 
right-handed zero modes for A; < 0. For k = there are 4 vector like zero modes for vanishing toron 
field, tfj^ = 0. For k > the zero modes have the form 

^f[x) = fix) toll" = ce^-('^(^')+^^(^)+*'^""/2)g{W2r)(.2_|.p)^(^^ ^e/o ^ 

with z = {xi + ix2)/Li — {t2 — iti)eL2/{2-Kk). H{z) denotes a polynomial in Jacobi's 0ri(-z) = &n{z\T)- 
functions, r = iL2/Li^ see Ref. 20. For k = n one has n linear independent polynomials H^^''^\z), 



see Table 2.4. The spin-flavor dependence of these zero modes is simply contained in the constant 



spin-flavor vectors They are chosen as eigenvectors of a chiral isospin transformation A (the 

main automorphism for differential forms in the geometric formulation), denoted by even and odd. 
With 7^ chosen diagonal one gets 

^Ib = ^a,l 4,1 , UJ°f, = 6a,l 6b,2 , Awli" = lla'^%ilh = ■ (A.5) 



For A; < one has to take the charge conjugated solutions Xab of Eq. ( A. 4 ). For A; = 0, = the 
zero modes still flt the form of Eq. ( |A.4| ) with H{z) = 1. However, the charge conjugated functions 
are then solutions at the same time, thus there are 4 zero modes in this case. 
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Table 2.4: The polynomials in ^-functions H^^'"'\i = 1, . . . ,n, for n = 1, . . . , 5. 



k = l: 






03 






k = 2: 




of 


01 






k = 3: 




of 63 


01 


Oi 02 Oi 




k = 4: 


Of 


el el 


01 


Ol O3 02 O4 




k = 5: 




el el 


01 


Oi ol O2 ^4 


of 02 O4 



The geometric version of these zero modes is given by 

X^'°{x,H)dx^ = f{x)oj%\ uj% = l-idx^^,oj"H = dx^ + idx^ . (A.6) 

In this form they can be easily compared with AZMs on the lattice. In the discretization of DK 
fermions to staggered fermions, the homogenous forms oc dx^ correspond to the staggered fermion 
field on the lattice sites with x^ = odd for ^ £ H. Since, the chiral isospin transformations generated 
by A in Eq. (|A.5D hold on the lattice, the lattice AZMs can be distinguished in even and odd 
modes too. The even AZMs then live only on the even lattice sites, xi + X2 = even. In regions of 
approximately constant f{x) they should be proportional to Q(x) 

^}{x) = for xi,X2 even , = a^^ = —ia? for xi,X2 odd . (A--7) 

Due to the particular spin-flavor dependence, the product rule of differentiation leads to a multi- 
plicative structure of the zero modes for A; > (A; < 0), proportional to the same spin-flavor vector oj. 
Let XAix) = fAix) a; be a zero mode in the background field A^{x). Then xa'Xb (x) := fA{x)fB{x) uj 
is a zero mode in the background Afj_{x) + -B^(x) 

r - ie{A^ + B^)] XA-XB = fs 7^ [9^ - ieA^] XA + /a 7^ [^^. - ieB^] Xb = ■ (A.8) 

For vanishing toron field = this multiplicative structure is evident from the explicit form of the 
zero modes, see Eq. ( A.4| ) and Table |2.4] . For other values of this structure is realized with the 



addition formulas for 0-functions.^'' 

A similar picture doesn't arise for a zero mode in the background A^{x) — B^{x), where A; > 
{k < 0) for both fields. Functions like give rise to a zero mode, if and only if they 

are continuous, which in general is not the case. Corresponding to that for k[A] > 0,k[B] < the 
corresponding spin-flavor vectors are complex conjugate to each other, and Eq. ( |A.§| ) is not valid. 

Finally we want to describe what happens, if a topologically non-trivial gauge field on T2 is 
embedded in a larger volume 7^'. Even if the gauge field on 7^' vanishes on most of this bigger torus, 
the zero mode does not become arbitrarily small for large distances to the region of non-vanishing 
gauge field. However, on the lattice we could numerically achieve a behavior as given in Eq. (|2.8|), 



see Figure 2.1. As a consequence of this non- localization of zero modes even in localized gauge field 
configurations, it is not clear how to construct (approximate) zero modes of instanton anti-instanton 
pairs. There is no simple way to glue the zero modes x X* of a local instanton (k = 1) and a 
local anti-instanton (A; = — 1) together, even if they are situated far from each other. In a transition 
region the spin-flavor given by lo must be fliped to uj*. It is not clear, whether such a spin-flip is 
possible without a significant rise of the eigenvalue. 
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